Statistical Estimation of Quantum Tomography Protocols Quality 
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A novel operational method for estimating the efficiency of quantum state tomography protocols 
is suggested. It is based on a-priori estimation of the quality of an arbitrary protocol by means 
of universal asymptotic fidelity distribution and condition number, which takes minimal value for 
better protocol. We prove the adequacy of the method both with numerical modeling and through 
the experimental realization of several practically important protocols of quantum state tomography. 
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Introduction. Quantum states and processes being 
a fundamental tool for basic research, also become a re- 
source for developing quantum technologies [l[: this de- 
mands for their characterization. At present quantum 
state/process tomography serves as a standard instru- 
ment [2| for characterizing quality of preparation and 
transformation of quantum states. Basically it includes 
a given set of unitary transformations over the state to 
be reconstructed, then the transformed state is measured 
and finally some computational procedure applied to the 
measuring outcomes completes the state reconstruction. 
The result is a state vector or a density matrix. In 
order to compare different schemes the achieved result 
can eventually be checked by calculating some distance 
measure, like fidelity, that shows the quality of prepara- 
tion/reconstruction of the quantum state. 
It is worth to mention that in real experiments the accu- 
racy of the reconstruction depends on two types of un- 
certainties: statistical and instrumental ones. If the total 
number of measurement outcomes (sample size or statis- 
tics) is large enough, the instrumental uncertainties dom- 
inate over fundamental statistical fluctuations caused by 
probabilistic nature of quantum phenomena [3|. Practi- 
cally the required statistics, allowing to exclude statis- 
tical fluctuations, depends on the tomographic protocol 
itself and the total accumulating time needed for taking 
data. From this point of view it is of the utmost inter- 
est to point out simple and universal algorithms for the 
estimation of the chosen protocol on design stage before 
doing experiments as well as the sample size for desirable 
quality of the state reconstruction. 

In this letter we propose a universal method for estimat- 
ing the quality of any tomographic protocol based on 
discrete degrees of freedom and test it with well-known 
reconstruction protocols of polarization states of qubit 
pairs. Since it represents a paradigmatic example, here 
we restrict ourselves to polarization degrees of freedom 
only. Manipulations with polarization qubits, qutrits and 
ququarts have been discussed widely in context of quan- 



tum state reconstruction and its optimization (3- 11|. 

Theory and numerical simulation An arbitrary s- 
dimensional quantum state is completely described by 
a state vector in a s-dimensional Hilbert space when it 
is a pure state, or by a density matrix p for a mixed 
one. To measure the quantum state one needs to perform 
a set of projective measurements. According to Bohr's 
complementarity principle, it is impossible to measure 
all projections simultaneously, operating with a single 
representative of the quantum state only. So, first of 
all, one needs to generate a set of copies of the state @. 
Then for each measurement j outcome rates Xj can be 
evaluated from amplitudes of a quantum process 



Xj = \Mi 



(1) 



The amplitudes Mj cannot be measured directly, but 
they are linearly related to the state vector For an 
arbitrary protocol based on m measurements the process 
amplitudes can be represented as: 



Mi 



X jC , j = l,2 t . 



(2) 



c being the state vector and Xj a row of the so called 
instrumental matrix X, which describes the entire set 
of mutually-complementary measurements. The normal- 
ization condition for the protocol bounds the total ex- 
pected number of events n and the acquisition time tf. 
jyjLi t r (Xj~Xjp)tj — n. We define the row of the mea- 
surement matrix B for a tomographic protocol as the 
direct product Bj = tjXj X* , its size being m x s 2 
and we assume m > s 2 . By using the matrix B, the 
protocol can be compactly written in the matrix form: 



T = Bp 



(3) 



Here p is the density matrix, given in the form of a col- 
umn (second column lies below the first, etc.). The vector 
T of length m records the total number of registered out- 
comes. The algorithm for solving equation (|3|) is based 
on the so called singular value decomposition (svd) fliij . 
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Svd serves as a base for solving inverse problem by means 
of pseudo- inverse (PI) or Moore-Penrose inverse [l4| • In 
summary, matrix B can be decomposed as: 



B = USV A 



(4) 



where U (to x to) and V (s 2 x s 2 ) are unitary matri- 
ces and S (to x s 2 ) is a diagonal, non- negative matrix, 
whose diagonal elements are " singular values" . Then © 
transforms to a simple diagonal form: 



Sf = Q 



(5) 



with a new variable / unitary related to p as / = V + p 
and a new column Q unitary related to the vector T as 
Q = U + T. We use this algorithm as a starting approx- 
imation for maximal likelihood state reconstruction. By 
defining q as the number of non-zero singular values of 
B we formulate two important conditions of any tomog- 
raphy protocol, namely its completeness and adequacy 
[15| . The protocol is supposed to be informationally com- 
plete if the number of tomographically complementary 
projection measurements is equal to the number of pa- 
rameters to be estimated; mathematically completeness 
means q = s 2 . Adequacy means that the statistical data 
directly correspond to the physical density matrix (which 
has to be normalized, Hermitian and positive). However, 
generally for mixed state it can be tested only if the pro- 
tocol consists of redundant measurements, i.e. to > q. 

PI provides with zero approximation for estimated pa- 
rameters of quantum states. Optimization of these pa- 
rameters can be done following the approach developed 
in Q in the frame of maximal likelihood method. It is 
worth to stress that the reconstructed state vector ex- 
tracted from the likelihood equation associates with fi- 
nite statistics of the registered outcomes of an exper- 
iment and therefore takes random values. The differ- 
ence between the reconstructed state vector and its ex- 
act value is caused by statistical fluctuations 11611 . The 
complete information matrix introduced in [3, [l2| allows 
analyzing arbitrary functions on fluctuating parameters 
of quantum states. One very relevant function in the 
context of quantum state/process estimation is fidelity 



F 



where po and p being the ex- 



act and the reconstructed density matrices correspond- 
ingly, tending to identity for complete protocols with an 
unlimited increase of a sample size. 

Basically the problem of the state reconstruction splits 
into two scenarios. The first one relates to the case when 
the reconstructed state is known in advance (like in quan- 
tum cryptography). Then statistical comparison between 
exact theoretical state and reconstructed one allows one 
either to reveal the source of instrumental uncertainties 
or to be sure that such uncertainties are small enough in 
comparison with statistical fluctuations. The second case 
occurs when the reconstructed state is unknown. Then 



using the approach developed below one should calculate 
the statistical distribution of fidelity between unknown 
exact state and reconstructed (by means of maximal like- 
lihood method) one. Fidelity will be either inside the un- 
certainty boundary for small statistics or outside for large 
statistics (when instrumental uncertainties prevail). In 
this case it is convenient to analyze the so called fidelity 
loss dF = 1 — F. When the reconstruction accuracy is 
determined by a finite number of representatives of the 
quantum state, fidelity loss is a random variable whose 
asymptotic distribution can be represented as [l7j |: 



(6) 



where £j ~ N (0, 1) are independent and normally dis- 
tributed random variables with zero mean values and unit 
dispersion, dj ~ — > 0, j max — 2s — 2 for pure states, 
jmax = s 2 — 1 for mixed states. The coefficients dj can 
be extracted from the complete information matrix [l7j |. 
Distribution ([6]) is a natural generalization of the x 2 - 
distribution for which all dj — 1. As it follows from © 
the average value of fidelity loss is: 



(1-F) = 



(7) 



PI implies introduction of so called condition number K , 
which is defined as the ratio between the minimal nonzero 
singular eigenvalue of B and maximal one 



K = 



(8) 



K determines the stability of the linear system Q and 
therefore can be used as a practical quantifier for esti- 
mating efficiency of the protocol: the lower K the better 
the protocol. If at least one of s 2 singular eigenvalues 
bj is close to zero then the protocol becomes incomplete 
and K — > oo. The optimal value of K would be unity 
that means uniform distribution of singular values [lij ]. 
The present approach allows to analyze arbitrary proto- 
cols of statistical reconstruction of quantum states. As an 
example we have selected three popular protocols. The 
first, suggested in Q and dubbed J16 Q, is suited for re- 
constructing 4-dimensional polarization states, as photon 
pairs degenerate in frequency generated in the process 
of spontaneous parametric down-conversion (SPDC). In 
this protocol the projective measurements upon some 
components of Stokes vector are performed for each qubit 
in pair individually, so the 16 two-qubit measurements 
are HH, HV, VV, VH, RH, RV, DV, DH, DR, DD, 
RD, HD, VD, VL, HL, RL where H, V, R, L, D de- 
notes horizontal, vertical, right and left circular and 45° 
diagonal polarizations respectively. Here, for example, 
the measurement setting HR means measuring horizon- 
tal polarization on the first qubit and right circular po- 
larization on the second qubit in pair. In our formalism 
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the protocol can be represented by an instrumental ma- 
trix with 16 rows Xj. 

In the case of independent measurement of two qubits the 
projective measurements can be chosen arbitrarily. In Q 
measured qubits were projected on the states possess- 
ing tetrahedral symmetry (R4) . There are several works 
showing that due to the high symmetry such protocol 
provides a better quality of reconstruction [H 0] . Let us 
stress that since m = s 2 = 16 the adequacy of both R16 
and J 16 can not be tested. 

We also consider a protocol where the whole single-beam 
two-qubit state is subjected to linear transformations us- 
ing two retardation plates [H (B144) . The total number of 
projective measurements is redundant and equals to 144, 
so the corresponding instrumental matrix has 144 rows 
and protocol admits adequacy testing, being m > 16. 
For these protocols we calculate condition numbers K 
which take the following values: Kris = 3, Kji 6 « 10, 
Kbih ~ 60[18|- Thus, we expect that the symmetrical 
protocol R16 provides with better state reconstruction 
quality [ij]. We have checked this statement with nu- 
merical simulations of each protocol applied to different 
two-qubit states depending on the sample size. 
As an example, let us consider the numerical reconstruc- 
tion of the Bell state |$") = -2= (|ififf 2 ) - \ViV 2 )). Fig. 
Q] shows average fidelities, calculated according to ([7]) 
as functions of sample size for each protocol. It turns 
out that the difference between curves disappears at suf- 
ficiently large sample size (10 5 ) 0, but for the same 
quality of state reconstruction the correct choice of pro- 
tocol allows using a smaller set of statistical data, i.e. 
finally reduces total acquisition time. Fig. Q] shows that 
protocols are ranged in accuracy as following: R16, J16, 
and B144 in complete agreement with the range given 
by condition number K. Figj2] presents accuracy distri- 
butions calculated according to ((6|) for the sample size 
3 x 10 3 . For better visualization we transformed the ab- 
scissa scale into particular common logarithmic one as 
z = —lg(l — F). Corresponding integers indicate the 
number of nines in fidelity recording. The density dis- 
tribution for the R16 protocol is narrower than the one 
for J16 and B144 and localizes in the region of lower 
losses or higher fidelities. The distribution for B144 is 
broader and lower in comparison with R16 and J16. Ob- 
viously, the narrower distribution of fidelity indicates a 
better reconstruction quality in the sense that the recon- 
struction procedure returns a better-defined state. Thus, 
FigE] confirms our expectation based on estimation of 
condition number K: R16 achieves the best results. 

Experiment. We prepared a family of two-photon po- 
larization states which can be easily converted into both 
entangled (in polarization) and factorized states 

|*) = (c 1 \H 1 H 2 ) +c 2 e llp \V 1 V 2 )) (9) 

with real amplitudes c\ and c 2 . The set-up is schemat- 
ically depicted in Fig. [3] A cw argon laser beam at 
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FIG. 1: (Color online) Dependence of the average fidelity on 
number of registered events forming the sample. 
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FIG. 2: (Color online) Density distribution of the scaled fi- 
delity z (lower abscissa) at n = 3 x 10 3 . Upper abscissa 
presents regular fidelity. 



A = 351nm pumps, after having selected by a Glan- 
Thompson prism (GP) its horizontal polarization, two 
type-I BBO crystals (1 mm) positioned with the planes 
that contain optical axes orthogonal to each other. The 
halfwave plate placed in front of crystals (A p /2) rotates 
the polarization of the pump by the angle (f> which con- 
trols real amplitudes c\ and c 2 in The crystals 
are cut for collinear frequency non-degenerate phase- 
matching around central wavelength 702nm. The rela- 
tive phase shift if in ([9]) is controlled by tilting quartz 
plates QP. If 4> = we prepare the state |^) = |Vi V2), 
if cf) = 22.5 and ip = 37r/2 then the state transforms to 
I*) = 71 d ff i ff a> - \ViV 2 )) = |$-). 

To maintain stable phase-matching conditions, BBO 
crystals and QP are placed in a closed box heated at fixed 
temperature. The lens L couples SPDC light into the 



4 




GP L/2 : QP BBO UVF WP WP, BS L GP M 



FIG. 3: (Color online) Experimental setup for different tomo- 
graphic reconstructions of photon pairs with variable polar- 
ization entanglement. 



monochromator M (with 1 nm resolution), set to trans- 
mit "idler" photons at 710nm. The conjugate "signal" 
wavelength 694 nm is selected automatically by means 
of coincidence scheme. Zero-order wave plates (A/2, A/4) 
are used in both J16 and R16 for arranging the projective 
measurement set. In B144 these plates were removed and 
additional achromatic quartz plates WPi^ (0.9183 mm, 
0.9167 mm) established the necessary measurement set 
at given wavelengths [j| . 

For protocols described above the statistical reconstruc- 
tion of prepared states ([9]) has been performed at given 
sample sizes. As an example, Fig.(Q} shows calculated 
widths of fidelity distributions at 1%- and 99%- quantiles 
for the three protocols (see Figj2]) as well as the experi- 
mentally reconstructed values for the specific state |$~). 
The approach described above provides the ideal accu- 
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FIG. 4: (Color online) Reconstruction of | $~^-state. Vertical 
bars show 1%- and 99% - quantiles for fidelity distributions. 
Dotted lines connected lower bar ends indicate critical signif- 
icant levels. 

racy level for quantum state reconstruction. It means 
that fluctuations of the estimated quantum states cer- 
tainly cannot lead to uncertainties smaller than this limit. 
Presence of instrumental uncertainties makes this level to 
be exceeded. Indeed, Fig|4] shows that, beginning from 
some sample size, the experimental value of fidelity falls 
out the theoretical uncertainty boundary shown as dotted 
lines (corresponding to 1% significance level which char- 



acterizes a given protocol). This happens since instru- 
mental uncertainties prevail over the statistical ones and 
indicates that either state preparation stage or measure- 
ment procedure were not performed accurately enough. 
Incidentally, the absence of uncertainty bars for data in 
Fig. 4 derives from running once the corresponding proto- 
col. Repeated measurements would provoke appearance 
of a statistical uncertainty, however respective theoreti- 
cal uncertainty levels should be re-scaled. 
Definitely comparison between reconstructed states and 
fundamental statistical fluctuations can serve to achieve 
precise adjustment of the set-up, detection of the unap- 
proved incursion into communication channel, etc. 

Conclusion. We have proposed and tested, both 
by numerical calculation and experimental realization, a 
new estimation scenario of tomographic protocols. Our 
results demonstrate the potentialities of this method for 
a widespread application to experiments on fundamen- 
tal quantum optics and quantum technologies. Indeed, 
it can provide in advance, based on condition number 
K , indications on the uncertainty that can be reached 
for a certain set-up by a tomographic scheme, providing 
experimentalist with a tool for choosing the best proto- 
col based on available experimental resources (retardant 
plates, polarization filters, etc.) and limited available 
time for data acquisition. Also we would like to stress 
that the developed method is quite general and that it 
can be applied to any sort of quantum states and mea- 
surement sets [15| . 
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